Static critical behavior in the inactive phase of the pair contact process 
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Steady state properties in the absorbing phase of the Id pair contact process (PCP) model 
are investigated. It is shown that, in typical absorbing states (reached by the system's dynamic 
rules), the density of isolated particles, pi, approaches a stationary value which depends on the 
annihilation probability (p); the deviation from its 'natural' value at criticality, p" at , follows a 
power law: p" at — pi ~ [p — Pc) 131 for p > p c . Monte Carlo simulations yield f3i = 0.81. A 
cluster approximation is developed for this model, qualitatively confirming the numerical results 
and predicting /3i = 1. The singular behavior of the isolated particles density in the inactive phase 
is explained using a phenomenological approach. 

PACS numbers: 05.70.Ln, 82.20.Mj, 64.40.Cn 



I. INTRODUCTION 



In the simplest models undergoing absorbing state 
phase transitions in the directed percolation (DP) univer- 
sality class, like the contact process (CP), the stationary 
state of the system in the inactive phase is the state de- 
void of particles. Other models in the same class, like the 
pair contact process (PCP), have however a richer struc- 
ture, associated with the existence of an infinite number 
of absorbing states [QJ^]. In the inactive phase, the ab- 
sorbing state the system evolves to depends on the ini- 
tial conditions and the distance from the critical point, 
and so does the average density of isolated particles in 
the stationary state. In the case of PCP, the field re- 
sponsible for the dynamics (the density of pairs of par- 
ticles, p 2 ) is coupled to another field (the density of iso- 
lated particles, p\). This background of isolated particles 
is responsible for the non-universality of some dynamic 
properties of the system at criticality |}|4|] . Recently the 
one-dimensional PCP with particle diffusion (known as 
PCPD or annihilation/fission model) has received a lot 
of attention and was at the center of some controversy 
H|]|||7|||; nevertheless its critical behaviour is not yet 
fully clarified. Carlon et al. || presented a pair approxi- 
mation study of this model which however is not valid for 
the non-diffusive case. The expected behaviour of PCP 
is qualitatively different from the D — limit of PCPD: 
in the non-diffusive case it is crucial to distinguish be- 
tween particles and isolated particles, since all activity 
stops as soon as the number of pairs vanishes, whereas 
isolated particles may produce a pair if particle diffusion 
is allowed. 

In PCP, as the system approaches the critical point 
from the active phase, the nonordering field approaches 
a 'natural' value (pi at ), and its behavior is described by 
the same power laws as those of the order parameter [ [To| . 

In the present work we investigate the behavior of this 
model in the inactive phase by looking at properties of its 
natural absorbing states, the absorbing states selected by 
the system's dynamics. As shown by both mean-field-likc 



approximations and Monte Carlo (MC) simulations, the 
stationary density of isolated particles develops a power 
law singularity as it approaches the critical point. 



II. PHENOMENOLOGICAL APPROACH 

In the inactive phase, the concentration of pairs is 
known to decay exponentially: /02(f) ~ e _Qt , with a ~ 
(p — Pc)" 11 • The exponent vu is associated with the tem- 
poral correlation length and mean-field theory predicts 

II = L . 

Now, if one is interested in studying the stationary sin- 
gle particle concentration as a function of p(> p c ), one 
can refer to the coupled Langevin equations describing 
this dynamical process jOL At the mean-field level, the 
time evolution of p\ is given by 



dpi 2 
—7- = np 2 - wipip2 - uip 2 
at 



(1) 



obviously coupled to the evolution of p 2 . With the change 
of variables pT = p\ 



equation (1) reads 



dn _ 2 

— = -W1P2P1 - uip 2 . 



(2) 



In the following, we will assume /02(f) 
Then, equation (2) can be solved exactly: 



P 2 (0)e- 



7T(*) = 
with 



"1" ( Wip 2 {0) e -crf + 1 

w\ \ a 



C^(0) + ^,2(0) + ^) 



In the limit f — ► 00, one has 

pi(t — > 00) = j- + G. 



"1P2(") 



(4) 



(5) 
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In the vicinity of the critical point, when a — > 0, C 
decreases more rapidly than a, and one can therefore 
conclude that 

pl(t->- oo) ~ (6) 
w i 

At the critical point, ~pi(t — > oo) = 0, i.e. 

ft ( t _>oo) = H ( 7 ) 

p™ a * is then given by r\/wi, as already found in [ p^[ . 
But now we can see that — > oo) varies with p, in the 
inactive phase: 

Pl (t ^ oo) ~ P r* ~ ^a. (8) 

We must point out that the present argument relies 
on a mean-field approach and on the assumption that 
P2(t) = P2(0)e~ at (which may not be true at early times). 
A variation of pi(t — > oo) with p is then plausible. How- 
ever, the prediction that the critical exponent of the 
quantity p\ at — p\ should be equal to v\\ is only valid 
within a mean-field approximation and is not expected 
to apply to the Id PCP model. 

III. CLUSTER APPROXIMATION 

Mean-field-like kinetic equations for the PCPD were 
obtained in || . Whereas the single-site approximation is 
only appropriate to the high diffusion limit (D — > 1), the 
pair approximation gives a good qualitative picture of 
the model for D > 0, but does not show some important 
characteristics of the PCP without diffusion. Indeed, ac- 
cording to the pair approximation, the steady state single 
particle density vanishes at the critical point (p c = 0.2 for 
the pair approximation with D = 0); however, it is well 
established that, at the critical point, the single particle 
density approaches a nonzero value, pi at (pi at — 0.242, 
in a sequential dynamical process). Also, the pair ap- 
proximation predicts a power-law temporal decay for the 
pair concentration in the inactive phase, contrary to the 
exponential decay found in the simulations. 

In the coarse-grained Langevin description, the fields 
that are used to characterize the system configurations 
are the local pair density and the local density of iso- 
lated particles. If one uses two-site clusters, as in the 
pair approximation, configurations with isolated parti- 
cles are not treated appropriately. Indeed, if one aims at 
building a better approximation, one has to go to three- 
site clusters and consider Poio, the probability of having 
an isolated particle at the center of the cluster. The 
price to pay is that of increasing complexity: the num- 
ber of variables and equations increases and these have 
to be solved numerically; Carlon et al. dealt with two 
variables and two equations, so they were able to obtain 
some results analytically. We have used a (3, 2) - cluster 



approximation |12| whose technical details are given in 
the Appendix. 

In Fig. I, we have plotted the stationary values of Pn 
and P\o as functions of p. As shown, the concentration 
of pairs vanishes for p > 0.128, i.e. p c — 0.128 within the 
present approximation. The estimate obtained by MC 
simulations is p c — 0.0771; so, as expected, the present 
result is an improvement when compared to what was ob- 
tained by the pair approximation. More interestingly, the 
concentration of particles at the critical point (pi = Pio , 
since Pn = 0) is nonzero; indeed P \ = 0.229, a value not 
far from p" at — 0.2418(2), as obtained by the simulations 
]l3"t . On the other hand, a linear decrease of p\ for p > p c 
can also be noticed. 



IV. NUMERICAL SIMULATIONS 



In this work, we concentrate on the critical behaviour 
of p\ in the inactive phase. The simulations that we 
present were done on systems of size L = 5000; to ensure 
a sufficient number of pairs we chose an initial parti- 
cle concentration of 0.5. The numbers of particles and 
of pairs were recorded versus time up to a number of 
Monte Carlo steps, t max , which ranged from t max = 10 6 
for p very close to the transition (p — 0.0775) down to 
tmax = 10 4 for p > 0.085. These times were chosen such 
that most of the samples (typically around 20000) had 
already entered the absorbing state. For each p we eval- 
uated the average activity time r av . To obtain the sta- 
tionary value of pi, we have averaged the final number 
of particles in those samples which entered the absorb- 
ing state at a time around T av . The short lived samples 
were excluded in order to eliminate finite size effects. On 
the other hand, in systems of this size, and for p close 
to p c , the concentration of pairs is still considerably high 
(typically around 0.1) at very long times - thus the con- 
centration of isolated particles is lower than in samples 
where the number of pairs is vanishingly small. For this 
reason, we have also ignored the contribution from these 
long-lived samples. We have chosen to consider samples 
with 0.75r av < r < 1.25T a „, but have checked that the 
results are stable with respect to other choices in a rea- 
sonable range. 

In Fig. 2 we show the stationary single particle con- 
centration as a function of p. Fig. 3 shows a log- log 
plot of (p^ at - pi) vs (p - p c ), with p\ at = 0.2418 and 
Pc = 0.0771. A linear fit of the data is clearly appropriate 
and leads to an exponent f3\ = 0.81(3). This seems to be 
a new exponent, not simply related to the DP exponents 
that describe the critical behaviour of other quantities, 
in the active phase of the PCP model. 
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V. CONCLUSION 

In this work, we show that the approach to the criti- 
cal point in the in inactive phase of the Id PCP model 
is signaled by a power-law singularity of a static quan- 
tity - the deviation of the particle concentration from its 
critical point value, p\ at . A cluster-mean- field approxi- 
mation was developed that predicts a finite value of Pi at , 
in reasonable agreement with MC simulations and a lin- 
ear decrease of p\ with (j> — p c ); this is confirmed by a 
phenomenological approach. Such a decrease in p\ is also 
shown in our MC simulations; the best fit to the data is 
consistent with an exponent (3\ = 0.81(3), different from 
the mean-field prediction. 

Critical behavior of a static quantity in the inactive 
phase has been observed by Lipowski and Droz in a 
rather different model with infinite number of absorbing 
states. In that case, a random walk argument leads to a 
relationship between the corresponding exponent and the 
order parameter exponent (3. In the case of PCP such an 
argument cannot be applied and whether /?i is related to 
the (DP) critical exponents that characterize the active 
phase of PCP requires further investigation. The same 
applies to other systems with infinitely many absorbing 
states. 
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3— site clusters in configurations 110,100 and 101 
are, respectively, Ay = +1, Az = +1, Av = 0. 

On the other hand, according to the PCP dynamic 
rules, a pair can also create a new particle at a ran- 
domly chosen nearest neighbour provided this is vacant. 
If one considers the configuration 1 1_J_ 11 (whose prob- 
ability is given, within the present approximation, by 

2 

An = X ^p ' w ith d = Pio = v + z), then, with rate 1 — p, 
the configuration 1 1 1 1 1 1 is generated. This change 
corresponds to Ax = +3, Ay = — 1, Az = 0, Aw = — 1 . 

The kinetic equations for x, y, z, v are obtained by con- 
sidering all the contributions from all the possible 6— site 
configurations: 

dx 

— = p[-4Ai - 6A 2 - 4A 3 - 2A 4 - AA 6 - 2A 7 - 2A 9 ] 

+ (1 - p)[3A n + Au + 3A 13 + 2A U + A 1B + A w 
+ 2A 17 +A 1S ) 

^ = p[Ai -A 4 -A 5 - 2A 7 - 2A 8 - 2A 9 - A w ] 
+ (l-p)[-A n -A 13 ] 

dz 

— = p[Ai + 2A 2 + 2A 3 + A 4 + A 5 + 2A 7 - A w ] 
+ {1 - p)[-A 12 - A 15 ] 

dv 

— = P [-2A 3 - 2A 5 - 2A 7 - 2A 8 ] 
dt 

+ (1 - p)[-A u + A u - A 13 - A u + A 15 - A 17 ) 
where 

A 3 2 2 2 

x x y x vy x y 

A\ = A 2 = — ^3 = „ ; A 4 = — ; 
c c J c d c 



APPENDIX y 2 v 2 x 2 yz xy 2 v y 2 zv 

A5 = 35-; Aq — 5—; A 7 — 5—; A 8 = 

cd^ c z d c z d cd z 

In the (3, 2)-cluster approximation |13| one uses 3— site 

clusters and allows for an overlap of 2 sites between ad- 2 2 2 2 

jacent clusters. Thus, the probability of a 6— site cluster ^ g = xy z ■ a w = V. z ■ ^ n = xy v ■ 

in the state ABCDEF is written as c 2 d ' cd 2 ' cd 2 ' 

□ _ d Pbcd Pcde Pdef , n s 

Pabcdef - Pabc ■ 0) x _yp y ^ 

Al2 C d(l-c-2dy 13 cd 2 ' 

We use x = P in , y = P no (= P011), z = P%oo(= P001), 
v = P101 as independent variables. It can be easily shown 

that P io = v + z - y, and P 00 o = 1 - x - y - 2v - 3z. _ xyvjv + z - y) _ ^ = y 2 z 2 

We then study all the processes that can occur when a 14 ~ cc p ' ~ cd(l — c — 2d) ' 

pair belonging to a 6— site cluster is selected. 

Take for example the configuration 111111. Ac- 
cording to the above approximation, the probability of ^ xyzil — c — 2d — z) y 2 v(v + z — y) 
this configuration is A\ — x 4 /c 3 , where c = Pn = x + y. 16 — cd(l — c — 2d) ' 17 _ cd 2 
When the central pair is selected, then, with rate p, the 
configuration 1 1 1 1 is generated. The variation in 

the number of 3— site clusters in the state 1 1 1 is then . y 2 z{\ — c — 2d — z) 

Ax — —4; analogously, the changes in the number of 18 — c d{l — c — 2d) 
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These equations were solved numerically. 
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FIG. 1. Pair density, p2 ( . 
approximation. 
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-), as a function of p, within the (3, 2)-cluster 
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FIG. 2. Stationary particle density (in the inactive phase), as a function of p. 
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FIG. 3. Log-log plot of (p" at - pi) vs (p — Pc), with p" at = 0.2418 and p c = 0.0771. The straight line is a least squares linear 
fit, with slopes 0.81. 
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